The influence of rotation in radiation driven wind from hot stars: 
New solutions and disk formation in Be stars. 
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3 ; ABSTRACT 

m 

The theory of radiation driven wind including stellar rotation is re-examined. 

After a suitable change of variables, a new equation for the mass loss rate is 

derived analytically. The solution of this equation remains within 1% confidence 

when compared with numerical solutions. Also, a non-linear equation for the 

\Q ■ position of the critical (singular) point is obtained. This equation shows the 

O 

point. For a stellar rotation velocity larger than ~ 0.7 - 0.8 Vb re aku P , there exists 

r£^ only one critical point, located away from the star's surface. Numerical solutions 

crossing through this new critical point, are attained. In these cases, the wind has 

a very low terminal velocity and therefore a higher density wind. Disk formation 
^— > ' 

in Be stars is discussed in the frame of this new line driven stellar wind solution. 
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1. Introduction 

The theory of radiation driven winds has been very successful describing the observed 
terminal velocities and mass loss rates from hot stars. After the pioneering work of Castor, 
Abbott and Klein (1975, hereafter CAK) who realized that the force due to line absorption 
in a rapidly expanding envelope can be calculated using the Sobolev approximation (Castor 
1974). They developed a simple parameterization of the line force and were able to construct 
an analytical wind model. 

Simultaneously and independently, Friend and Abbott (1986) and Pauldrach, Puis and 
Kudritzki (1986) (hereafter FA and PPK, respectively) calculated the influence of the finite 
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cone angle correction on the dynamics of the wind. They found a better agreement between 
the improved or modified theory (hereafter m-CAK) and the observations for both M (mass 
loss rate) and v^ (terminal velocity) in a large domain in the Hertzsprung-Russell diagram. 
Furthermore, Kudritzki et al. (1989, hereafter KPPA) developed analytical formulae for the 
localization of the critical point, mass loss rate and terminal velocity with an agreement 
within 5% for v<x> and 10% for M, when compared with numerical calculations. 

The influence of the star's rotation was investigated by Castor (1979) and Marlborough 
and Zamir (1984). Both studies concluded that the effect of adding a centrifugal force term 
in the CAK equation results in a lower terminal velocity wind, but the mass loss rate is not 
substantially affected. Similar results were founded by FA and PPK. 

Concerning Be Stars, the m-CAK theory gives a good description for the polar wind, 
but it fails to account for a slowly accelerating flow, like the one in the equator of Be stars. 
Therefore several additional mechanisms has been proposed, based on the radiatively driven 
wind theory to explain the equatorial flow of these objects. The influence of magnetic fields 
and rotation were studied by Friend and McGregor (1984, see review from Cassinelli 1998). 
A model with azimuthal symmetry, rotation and viscous force was developed by de Araujo 
et al. (1994) and Koninx & Hearn (1992) incorporated sound waves in the wind dynamics 
of Be stars. A model driven mainly by thin lines has been proposed by de Araujo (1995), 
he obtains an outflow with a shallow expansion and a large mass flux when the line force 
parameters are considered as free and not in a self-consistently mode (Pauldrach 1989). A 
2-D hydrodynamical model from Bjorkman and Cassinelli (1993) and Owocki et al. (1994) 
concluded that a meridional current may be responsible for the concentration of matter 
towards the equator, forming a wind compressed disk ("WCD"). An assumption of the 
WCD model, that the line-force is strictly central was released by Owocki, Cranmer & 
Gayley (1996). They conclude that non-radial line-force components together with gravity 
darkening effects can inhibit the formation of a WCD structure. 

Despite all the efforts, the most important problem that pervades all these models, is 
the strong equatorial expansion that they exhibit. Not only the terminal velocity is too high 
(~ 1000 km sec" 1 ) but also there is a sharply increase of the velocity field. The fitting of 
observed H a line profiles requires terminal velocities of about 200 km sec~ l or less (Poeckert 
& Marlborough 1978), while from Fe II line profiles, Hanuschik (1994, see also Waters & 
Malborough 1994) concluded that the disk expansion velocity must not be much larger than 
the Doppler width. 

When rapid rotators are theoretically studied (such as Be stars), many authors have 
reported the appearance of numerical problems when the rotational speed is about 0.8 times 
the break-up speed, see e.g., FA, Poe and Friend (1986), de Araujo and de Freitas Pacheco 
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(1989), Boyd & Marlborough (1991). 

In view of the previous comments our purpose is to perform a re-analysis including the 
rotational centrifugal force term in radiation driven wind theory. 

In section 2, an analytical treatment with the inclusion of the star's rotational speed is 
carried on. In section 3, after a suitable coordinate transformation, exact formulae for the 
location of the critical (singular) point (s) and for mass loss rate are obtained. The roots of 
the critical point function defines the number of singular (critical) points and their location. 
We show the existence of a new family of singular points in addition to the standard one (m- 
CAK solution family). A simple, approximative, treatment for the location of the singular 
point (s) and the value of the corresponding mass loss rate (eigenvalue) is introduced in 
section 4. Furthermore, in this section, numerical and analytical results for the standard m- 
CAK solution are compared, showing the high confidence of this analytical approximation. 
Section 5 is devoted to find a numerical solution from the momentum equation of the wind, 
that starts at the stellar and surface reaches infinity after passing through a critical point 
that belongs to this new family of singular points. Also a comparison between numerical 
and analytical result are performed in order to show the accuracy of this analytic treatment. 
Finally the applicability of this model to explain disk-formation in Be-Stars is discussed in 
section 7. Conclusions and future lines of research are summarized in section 8. 



2. Hydrodynamic formulation 

The standard model for radiation driven stellar winds considers one component isother- 
mal fluid in a stationary regime with spherical symmetry, neglecting the effect of viscosity, 
heat conduction and magnetic field (CAK, FA, PPK). The continuity equation reads: 

47vr 2 pv = M (1) 

and the momentum equation is given by: 

dv I dp GM(l-T) uj(r) 
v— = — K — '- + -? — + g Hne (p, dv/dr, n E ), (2) 

dr p dr r z r 

here v is the fluid velocity, dv/dr the velocity gradient, p the mass density, M is the star's 
mass loss rate, p the fluid pressure, v^ = v rot R*/r, where v rot is the star's rotational speed 
at the equator, T is the radiative acceleration caused by Thomson scattering in terms of 
gravitational acceleration and g hne (p, dv/dr, he) is the acceleration due to the lines. The 
standard parameterization of the line force (Abbott 1982, PPK, FA) is: 

line C nv!„. .. A„IA„\ ( Jl„ dl! V ( H E ' " 



g^ = - 2 CF(r,v,dv/dr)^v-) [^ ) ■ (3) 
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The coefficient C is given by: 



C = TGMk [ ^-^- ) , (4) 

K a E v th M, 



here Vth is the thermal velocity of the protons, a E is the Thomson scattering absorption 
coefficient per density and n E is the electron number density in units of 10 _11 cm~ 3 , 



W(r) = -{1 - y/1 - (Rjr) 2 } (5) 

is the dilution factor, CF is the correction factor (see appendix B) and all the other quantities 
have their usual meaning (see, e.g., PPK). 

Introducing the following change of variables 

u = —R*/r, (6a) 

w = v/a, (6b) 

w' = dw/du, (6c) 

where R* is the star's radius and a is the isothermal sound speed, i.e., p = a 2 p. 
The momentum equation (2) with the line force (3) becomes: 



i) W ^ + A + I + ^ ~ C ° F 9{U){W) ~ 5 { W ^y 



F{u,w,w')= [1^ — )w— + A + - + a 2 rot u-C CFg(u){w)- d [w—) =0 (7) 
here 



_ GM(1 - p _ vj c 
A ~ a 2 R* ~ 2a 2 ' (8&) 



6 



and 

a ro t = — , (8d) 

a 



where v esc is the escape velocity, and D is defined by: 



(i + z ge y ge ) i 

(i + A He r He ) m H ' { ) 
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here Yn e is the helium abundance relative to the hydrogen, Z# e is the amount of free electrons 
provided by helium, Ajj e is the atomic mass number of helium and m# is the mass of the 
proton. 

The standard method for solving this non-linear differential equation (7) together with 
the constant C (the eigenvalue of the problem) is imposing that the solution pass through 
the singular (critical) point (see CAK for details), together with a constrain on the optical 
depth integral: 

/ a E p{r)dr = - (10) 

JR, 6 

Other authors prefer to use an equivalent lower boundary condition, e.g., setting the 
the density at the stellar surface to a specific value, 

p(R*) = p*, (11) 

de Araujo et al. (1989) use p* = 10~ n g cm~ 3 , Friend & MacGregor (1984) use p* = 1.5 x 
10~ 9 g cm~ 3 . 

A critical point is located where the singularity condition, is satisfied: 

^Lf(u,w,w')=0 (12) 

At this specific point, regularity is imposed, namely: 

d w A dF dF . 

- Fiu , w , w i ) = - + - W . =0 (13) 



3. Location of the Critical Point(s) 

In this section, after a coordinate transformation we obtain a set of analytical equations 
for the eigenvalue and the location of the singular point. For the last one, we show, that 
a new family of singular points exist. We also prove, that the standard m-CAK solution 
vanishes for stars with high rotational velocity. 



3.1. Coordinate transformation 

In order to solve equations (7), (12) and (13), we use the following coordinate transfor- 
mation: 

Y = w w' (14a) 
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Z = w/w' (14b) 

with these two new coordinates, the equations turns to: 



l ~Yzj Y +A + 2/u + a 2 rot u -C'f 1 (u,Z)g(u)Z' 5 / 2 Y a ' 5 / 2 =0, (15a) 

{ 1 -Yz) Y - C 'h{u, Z) g{u) Z- & ' 2 Y a - 5 / 2 = , (15b) 

{ 1 + Y^) Y -ZZ/u 2 + al ot Z -C'Uu,Z)g(u)Z- s / 2 Y a - s / 2 =0, (15c) 

derivation details and definitions of f\(u,Z), f 2 (u,Z) and f 3 (u,Z) are summarized in ap- 
pendix A. 

Solving for Y and C" from the set of equations (15a), (15b) and (15c), we obtain: 

1 , ( h \ ( , , 2 . 2 



and 

C'(M) = A-(i- —) z 5 ' 2 Y l - a+5 ' 2 (17) 

9J2 V YZ J 

These equations are generalizations of the relations founded by KPPA (see their eq. [21] and 
[34] for Y and eq. [20] and [44] for the eigenvalue) including now the rotational speed of the 

star. 



3.2. The critical point function R(u, Z) 

It is not possible to obtain the location of the critical point from this set of equations be- 
cause we have only three equations and four unknowns (Y,C, Z and u). But from equations 
(15a), (15b) and (15c), we obtain a function, R(u,Z), defined by: 

2 1Z / 2 \ 

R(u, Z) = -- + — - a 2 ot Z + f 123 (u, Z)(A + - + a 2 rot u (18) 



where f 123 (w, Z) is defined by 

f ( va _ f 2 (u, Z) - f 3 (u, Z) 

f 123 (u, Z) = — — — — — — 19 

h{u,Z) - Ji{u,Z) 

The root(s) of this function R(u, Z) gives the location of the critical (singular) point(s) u cr i t . 
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Notice that no approximation whatsoever has been used in the derivation of the above 
equations. 

In order to know the range of the variable Z for a typical hot-star wind, we have per- 
formed a full numerical calculation, i.e., we solve without any approximation: the momentum 
equation (7) together with the singularity condition, equation (12), the regularity condition, 
equation (13), and a lower boundary condition: equation (10) or equation (11). Our typical 
star is an 05 V star with the following stellar parameters: T e ff = 45000i^, log g = 4.0, 
R/R & = 12 and v ro t = 0, while the line force parameters are: k = 0.124, a = 0.64 and 
5 = 0.07. Figure 1 shows the behavior of Z, which ranges between and 2, for the whole 
wind. This figure also shows two /3-field approximations (see below). 

Knowing now the range of the variable Z we may plot R(u, Z) in terms of u and Z. 
Figures 2 and 3 show two different visualization of the function R(u, Z) for this 05 V star 
with v ro t/vbreaku P = 0.5. Figure 4 shows, besides R(u, Z), the zero-plane defined by R = 0. 
Thus, the root(s) of R(u,Z) is(are) the curve(s) defined by the intersection of this surface 
R(u, Z) with the zero-plane. Figure 4 shows both surfaces, their intersections are two families 
of critical points. The standard m-CAK family of solutions (CAK: locus of singular points) 
is located in the zone defined approximately by: u G [— 1, —0.3] and Z G [0, 3]. 

When full numerical calculations are carried on, the lower boundary condition, eq. (10) 
or eq. (11), fixes one point in this family of critical points and in this way an unique numerical 
solution for the m-CAK wind is achieved. 

As the rotational velocity of the star increases, this m-CAK locus does not longer inter- 
sects the zero-plane because the function R(u, Z) becomes negative in this region, as figure 
5 shows for the case v rot /vb rea k up = 0.9. Thus, standard m-CAK solution does not exist for 
stars with highly rotational speeds, e.g., Be-Stars. This is the reason why many authors 
(see .e.g., de Araujo et al. 1994, Boyd & Marlborough 1991, FA) have reported numerical 
difficulties in finding the location of the singular point for high rotational velocities, it does 
not exist. 

Furthermore, inspecting carefully figure 5, a new family of solution is found in the 
region defined approximately by u G [—0.2,0] and Z G [0,3]. This family of solutions has 
been always present, even for lower rotational velocities, see the case of figure 4. 

When both families of solutions are present, e.,g., for not highly rotational speeds, the 
family that satisfies the lower boundary condition, eq. (10), is the m-CAK locus. A study 
about the existence of a physical solution, that crosses through a singular point that belongs 
to the new family, when both families are present, is the scope of a forthcoming article. 



We are now ready to look for numerical solutions, that starts at the stellar surface, then 
crosses through a singular point located at this new locus and reaches infinity, for the special 
case when only the new family is present, e.g., for highly rotational speeds. We will perform 
this task in section 6. 

Before, in the next section, a simple approximation for the function R(u, Z) is intro- 
duced. This approximation allows an easy location of both, the standard and the new 
singular point. 



4. Approximative Treatment 

In this section we develop an approximative treatment of the critical point function 
R(u,Z). After introducing a /3-field approximation, the function R(u,Z) transform to 
R app (u). The zeros (roots) of this approximative function are the location of the critical 
points. Knowing the location of the critical point (s) it is straightforward to calculate the 
eigenvalues or equivalently the mass loss rate of the star. This approximative calculations 
are in much better agreement with full numerical calculations than previous analytical treat- 
ments and includes the influence of star's rotation. 



4.1. The approximative critical point function R app (u) 

From numerical calculations, PPK found that the velocity field, for stars with effective 
temperatures between 40.0007^ and 50.0007^ can be approximated to the so called /3-field 
approximation, namely: 

v = v OQ (l + uf (20) 

with (3 = 0.8. This relationship is broadly used for stellar wind diagnostics and it is justified 
a posteriori by the quality of the results achieved (Kudritzki & Puis, 2000). 



Applying this approximation for the Z variable, it becomes: 

„ (1 + u) 



P 



(21) 



Figure 1 shows the Z vs. u profile from numerical calculations and the behavior of eq. (21) 
for two different values of the f3 parameter ((3 = 0.8 and f3 = 1.0). This approximation holds, 
for the Z variable, up to ~ 2.5/?* (u ~ —0.4) for j3 = 0.8 and up to ~ l.li?* (u ~ —0.8) for 
P= 1.0. 
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Thus, replacing eq. (21) in eq. (18), R(u,Z) transform to R app (u), defined as: 
213 2(1 + u) o (1 + u) , , , / „ 2 „ \ 

Ra PP (u) = ~JY^ + U 2p ~ a rot-Jf^ + /l2 3 («) \A + - + U^uj (22) 

The function f\23(u,Z) transforms to fi23(u), and its behavior is shown in figure 6 for 
different values of the free-parameter (3. Once the root(s) of R app (u) in the interest domain 
(interval — 1 < u < or R* < r < oo) is(are) obtained, the values of Y and C = C'(M) are 
obtained from equations (16) and (17) respectively. 

4.2. An Example: CAK with Rotation 

It could be thought, that this simple case has only an academic interest, but as we shall 
see later, that these results contribute to the understanding of the dynamic of rapid rotators. 

For the CAK model (with 5 = 0), the f's functions are: f\ = 1, f'2 = f'3 = a and 
f 123(11) = 0, then the R app {u)(= Rcak(u)) function reads: 

Rcak(u) = - irT ^ y + ^-^ - a rot -^ (23) 

this is a fourth order polynomial in u. Contrarily to the expected, the non-rotational case, 
Rcak{u) does not depend neither on the stellar parameters: T e ff, M*, i?*, nor on the line- 
force parameters (k, a, 5), it explicitly depends on (3. Furthermore, for the rotational case 
it depends on the stellar parameters and the free-parameter j3, but not on the line-force 
parameters. 

CAK and KPPA have shown that the behavior of the velocity profile is given approxi- 
mately by: 

v(u) = VooVTTu (24) 

This results was obtained from approximate solutions for a non-rotating CAK wind, 
(see CAK eq.[47] or KPPA eq.[24]). This relationship corresponds to a /?— field with j3 = 
1/2. Figure 7 shows the function Rcak(u), for this 05 V test star, for different values 
of v ro t/vbreaku P and (3. From this figure, it is clear that only one critical point exists in the 
domain of interest. Bjorkman (1995) arrived to the same conclusion in his topological study 
of the non-rotating CAK model. 

As the rotational speed increases, figure 7 shows that the location of the CAK singu- 
lar point is shifted downstream, this is the same conclusion reached by Castor (1979) and 
Marlborough & Zamir (1984). 
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The zero {u cr u) of equation (23) in the interval — 1 < u < is giving by the analytical 
formula: 

Ucrit = -{h - h - I) (25) 



where the t-coefficients are: 



to = (-2 + a 2 rot + 2[3 2 ) (26a) 

ti = A32a 2 rot -2l6a 4 rot + 2l6a 2 rot t + 2t 3 (26b) 



3 



«i + V- 4 *o + *i J ( 26c ) 

2h 2 



3 a rot t 2 



3 a rot 2 3 



(26d) 

(26e) 



-8 + ™- + *£ (26f) 



°Vot °rot 



_2_*o 

3 a 



*6 = t^T" ( 26 ^ 



rot 



t 7 = VI - *e + *3 + U (26h) 

t 8 = j2-2t 6 -t 3 -U + -^- (26i) 

For the simple case of a ro t = 0, a straightforward and exact solution for equation (18) is 

u c = -Z c (27) 

and after applying the approximation Z = (1 + u)//3 with (3 = 1/2, the critical point 
is founded at u c = —2/3 or r c = 1.5 R*, as CAK and KPPA previously obtained using 
different approaches. 

Table 1 shows the accuracy of this analytical approximation. Columns 2-4 and 6-8 
show, for different values of (3 parameter, the location of the singular point (in terms of i?*), 
and mass loss rates (in units of 10~ 6 M & year~ l ) respectively. For comparison, columns 5 and 
10 give the values of r crit and M from full numerical calculations 1 . For the non-rotational 
case, the best values is j3 = 1/2 confirming the results of previous analyses. KPPA cooking 
recipe for Mqak (see KPPA eq. [27]) gives Mqak = 3.074 1O~ 6 M year^ 1 , same as the non- 
rotational case with /3 = 1/2. While the value of the mass loss rate is almost independent 



1 Here, full numerical calculation solves equations (7), (12), (13) and the lower boundary condition 
equation (10). 
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of the value of (3, the location of the critical point is very sensitive to the value of (3, but in 
the presence of rotation, this sensitivity disappears. The reason is that here the rotational 
term (third term in equation 23) is the dominant one, and therefore the dependence on (3 
becomes minimal. Furthermore, the larger the rotational speed, the larger is the (3 value 
that best fits the numerical calculations. 

The terminal velocity of the wind, v^, can be now computed from direct integration of 
equation (16). This follows because Y is the unique solution of (7) or (15a). 



4.3. The m-CAK Standard Solution 

Now we analyse the influence of the finite cone- angle effect. The main difference between 
this approach and the one of KPPA is that they introduced the /?— field approximation before 
calculating the derivatives of the singularity and the regularity conditions. Therefore they 
could not obtain the equations (15a), (15b) and (15c). In fact, the last two of them become 
the same in the KPPA approximation. 



4.3.1. OB stars 

In order to obtain the critical point and the mass loss rate it is necessary to know the 
value of (3. As discussed above, PPK found that a value of (3 — 0.8 is appropriate for OB 
Stars. In figure 8, same as fig. 1 but here in the region near the photosphere, the '+' symbol 
indicates the location of the critical point at the Z(u) profile. The /?— field approximation 
with (3=1 intersects almost exactly this location showing that for this type of stars (3=1 
is a fairly good approximation to calculate the critical point location and the mass loss rate. 

Table 2 shows a comparison between analytical and full numerical calculations for the 
05 V test star (see stellar and line-force parameters in section 3.2). The best agreement for 
both analytical, r cr u and M, with the numerical calculation occurs for (3 = 1. Furthermore, 
this predicted analytical values agree much better with the numerical calculations than the 
values given by the 'cooking recipe' from KPPA. This recipe gives for the non-rotational 
case: r crit = 1.038 R* and M = 1.976 10~ 6 M Q year' 1 (with (3 = 1). When the analytic 
approximation of KPPA is compared with numerical calculations, it is found that the values 
of M are always smaller than the numerical values. For the non-rotational case, our approx- 
imation match almost exactly the values of the numerical calculations, but for rotational 
cases it still shows values slightly less than the numerical ones. 
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4-3.2. Cental Stars of Planetary Nebulae 

Radiative driven stellar wind are also present in Cental Stars of Planetary Nebulae 
(CSPN), i.e., for objects where the effect of photospheric extension is important. Pauldrach 
et al. (1988) have computed radiative driven wind models along evolutionary tracks for these 
objects. They found that the the predicted terminal wind velocities are in agreement with 
the observational data, while the mass loss rate are in qualitative agreement. 

In their non-LTE analysis of CSPNs, Kudritzki et al. (1997) used a /?— field approxi- 
mation for the velocity, finding that (3 — 1.5 is the best value for these stars (see table 1 
from Kudritzki et al. 1997). Table 3 shows a comparison between analytical and numerical 
calculations for one of the models of Pauldrach et al. (1988). The stellar parameters for 
this model are: T eff = 80000K, logg = 5.24, M/M & = 0.565 and R/R Q = 0.3, line-force 
parameters are: k = 0.053, a = 0.709 and 5 = 0.052. This analysis confirms that a value of 
j3 = 1.5 is a fairly good one for non-rotational CSPN, but as the rotational velocity increases 
a slightly larger value must be assumed. 

In this section, we have developed an analytical treatment for the location of the criti- 
cal point (standard m-CAK) and mass loss rate that improves (includes rotation) previous 
investigations, and has about 1% confidence when is compared with full numerical calcula- 
tions. 



5. Analytical Treatment of the New Wind Solution. 

In previous sections, the existence of a new family of solutions, that are different of the 
standard m-CAK, has been established from the properties of the R(u, Z) function. In this 
section we study the approximate function R app (u) to get information about the location of 
the singular point and the mass loss rate. 



5.1. Critical Points 

5.1.1. The Standard m-CAK Critical Point 

Figure 9 shows R app (u), with (5 — 1, for an 05 V star for different values of the star's 
rotational speed for both models, CAK and m-CAK. From this figure it is clear that the 
inclusion of the finite disk correction factor in the wind momentum equation (7), 'twist' the 
Rcak( u ) function to R app (u). It is evident from this result that while Rcak(u) displays only 
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one root in the interest domain, the number of roots of R app {u) depends on the rotational 
velocity. 

For non-rotational cases (figure 9a) the CAK critical point shifts to a location close 
to the stellar surface, this is the standard m-CAK critical point. As the rotational speed 
increases, the location of the m-CAK singular point is shifted downstream (FA, PPK), but 
differently as in the CAK case, it occurs up to certain rotational speed. 

In order to have a better picture of the approximation Z — (1 + u)/f3 we have to recall 
the function R(u, Z) (see, e.g., figure 4). This is a surface in the phase-space defined by the 
independent coordinates u and Z. The approximation Z = (1 + u)//3 is a vertical-plane in 
this phase-space that cuts the i?-surface defined by equation (18). The projection of this 
intersection-curve over the vertical plane u — R is the R app function showed in figure 9. 

Figure 10 is a contour plot of the i?-surface at the zero-plane, the u — Z plane. Both 
continuous-lines correspond to the two families (loci) of singular points, discussed in section 
(5.1), while the dashed-line is the /?— field, approximation (/3 — 1). This approximation cuts 
the standard m-CAK family in two points. The standard m-CAK solution passes through 
the first intersection point (nearest to stellar surface, u = —1) and not through the second 
intersection point, as a result of imposing the lower boundary condition equation (10). The 
second intersection point is only a consequence of the /?— field approximation and a solution (if 
exists) that passes only through this point probably have no physical meaning for radiation 
driven winds. 



5.1.2. The New m-CAK Critical Point 

In addition to the standard m-CAK critical point, figure (9) and figure 10 shows a 
second critical point (root of R app ), the last intersection-point. As the rotational speed 
increases, for this 05 V star, the standard m-CAK critical point now disappears of the 
integration domain when v rot /vb r eaku P ^ 0.9 (fig. 9d, see also figure 5). In this case, the 
CAK critical point is shifted far downstream in the wind. At the critical point, u cr u, the 
function /123 {u cr u) — > 0, so R app (u) is almost the same for CAK and m-CAK models (see 
equations 22 and 23). Thus, for the fast-rotational case, equation (25) gives a very good 
approximation for the location of the critical point and therefore the mass loss rate, too. 

We conclude that the CAK model with rotation is applicable for high-rotational veloc- 
ities to calculate the mass loss rate and the location of the singular point of the m-CAK 
model. 
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Because we are now interested in stars with high rotational velocities, from now on our 
test star will be a typical B\ V Star with following stellar parameters: T e jj = 25000/^, 
logg = 4.03, R/R Q = 5.3 (Sletebak et al. 1980) and line force parameters: k = 0.3, a = 0.5, 
5 = 0.07. For this star, figure 11 shows R app (u) for two different values of f3. In both cases, 
as the rotational speed increases, R app (u) goes to a scenery where there is only one critical 
point, the new one. For (5 = 2.5, a lower value of the rotational speed is required to get to 
this situation, i.e., v rot /v breakup > 0.6, while for (3 = 1, the required is v rot /v breakup > 0.8. 
This last value is lower than the required for the 05 V Star (see fig. 9). 

In order to know all the families of solutions, a topological analysis is needed, similar 
to the analysis performed for the non-rotating CAK model by Bjorkman (1995). As a first 
step in order to learn about this new solution, we carry on numerical calculations only for 
the fast-rotational case, i.e., the case where only one root of the R app function exists, the 
new critical point. It is not clear yet, that a wind solution should cross only through one 
critical point or there are solutions that have to cross through more than one critical points 
as, e.g., in the magnetic wind model of Friend & MacGregor (1984). 



6. Numerical Results 

It has been established in the previous sections that in a fast-rotational case, the stan- 
dard m-CAK critical point vanishes and there exist a second critical point. In this section 
we will obtain full numerical solutions, that starting at the star's surface reaches infinity 
passing through this new critical point. 



6.1. The Numerical Method 

In order to solve the non-linear wind momentum differential equation (7) for a fast- 
rotational case, we will not use the 'standard' numerical method, i.e., once the singular point 
location, its velocity and velocity gradient are approximatively obtained, direct integration 
up and downstream are performed (e.g., using Runge-Kutta or Burlish-Stoer methods). 
After this process has been done, the lower boundary condition is calculated and the whole 
procedure starts again until convergence is achieved. 

Instead, we use here a finite-difference method, modified for handling singular points 
(Nobili & Turolla 1989). This method uses a trial solution as its initial guess. This initial 
solution then relaxes, using the Newton method, towards the numerical solution (see Nobili 
& Turolla 1989 for details). The main advantage of this method is the fact that it is not 
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necessary to give a priori the location of the critical point, as the standard method requires, 
so no guess value for (3 is needed for solve the momentum equation. 



6.2. The Fast Solution 

First we calculate the non-rotational case, v rot = 0, for our test star using r = 2/3 (eq. 
10) as lower boundary condition. In figure 12 the velocity profile for this cases (dashed-line) 
is shown. Numerical results are: critical point location at r cr u = 1.014 i?*; mass loss rate 
M = 3.178 10~ 9 M Q year" 1 ; and terminal velocity t>oo = 2025 km sec -1 . 



6.3. The Slow Solution 

We have solved numerically the fast-rotational case for v rot /vb rea k up = 0.8. This new 
solution is denser than the fast solution (m-CAK) and we call it hereafter slow solution 
due to the lower value of the terminal velocity achieved, figure 12 also shows this profile 
(continuous-line). The location of the critical point is at r crit = 23.63 /?*, confirming the 
approximative calculations of the preceding sections, that the critical point is located far 
away in the wind. The terminal velocity of this new solution is i^ = 443 km sec -1 . This 
value is a factor 4 lower than the standard m-CAK solution (the fast-solution, dashed-line 
in figure 12). The mass loss rate is M = 6.854 1O~ 9 M year -1 , approximately twice the 
value of the non-rotational case. 



6.3.1. The Approximation Z = (1 + u)/f3 

Figure 13 shows Z versus u from full numerical calculations (continuous line). In ad- 
dition, the location of the critical point is drawn with a '+' symbol in this figure. The 
approximation Z — (l + u)//3 (a straight line in this figure) is also plotted for two different 
values of the parameter (3. The line with (3 = 2 intersects Z(u) close to the location of the 
critical point, showing that for fast-rotational cases a value of (3 > 1 matches better the 
location of the critical point and the mass loss rate (eigenvalue). 
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6.3.2. Functions R(u, Z) and R app {u) 

In order to validate the R app {u) function, we have plotted in figure 14 the function 
R(u, Z) from the results of the full numerical calculations (continuous-line) and compared 
it with the R app (u) functions for two different values of the f3 parameter (dashed-lines). 
The behavior of R(u, Z) and R app {u) in the region far away in the wind (u > —0.05) is 
similar for both values of j3, confirming that, also for this new solution, the last root of the 
approximate function R app (u) gives an accurate value for the location of the critical point, 
r cr it, and through eq. (16) and eq.(17) the eigenvalue, C'(M), for the fast-rotational case. 



6.4. Vqo and M from slow solutions. 

As de Araujo (1994) and Stee at al. (1995, and references therein) argued, it is possible to 
obtain equatorial solutions with lower terminal velocities when the wind is driven mainly by 
thin lines, e.g., using a lower value of the line force parameter a. Thus, we have investigated 
the value of v rot that makes that the standard m-CAK critical point disappears remaining 
only the slow solution. 

Table (4) summarises the value of this minimum rotational velocity from the R app 
function for the existence of the slow solution. The values of a and k are taken from de 
Araujo(1994) and 5 = has been used. The values of v^ and M are coming from full 
numerical solutions. This result suggests that in order to obtain a better modeling of this 
slow-solution a self-consistent calculation (Pauldrach 1989) of the line force parameters is 
necessary. This table (4) shows, that the thinner the lines drive the wind (lower a), the lower 
is the required rotational velocity for slow solution and therefore the lower is the terminal 
velocity and mass loss rate. This table also shows a weak dependence of the /3-parameter 
on the minimal rotational velocity. 



7. Disk Formation in Be Stars. 

For our model is natural to think that a fast-rotating star can have both wind solutions, 
e.g., Be Stars. The polar direction is equivalent to the non-rotational case, while the equa- 
torial flow corresponds to the fast-rotational case. In figure 15 the density profile for both 
directions, i.e., polar (dashed-line) and equatorial (continuous-line) are plotted against the 
inverse radial coordinate u. 

The ratio between the equatorial and polar densities are about a factor 100 very close to 
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the star's surface but mainly a factor 10 for the whole wind. This density quotient value is 
still a factor 10 smaller than the observational values. The disk aperture angle occurs when 
the m-CAK critical point disappears from the R(u, Z) function, so this angle will depend 
strongly on the assumed value of the line-force parameter a. A detailed '2-D' (latitude 
dependent) velocity and density profile and line formation is the subject of a forthcoming 
article. 



8. Conclusions 

After a suitable change of coordinates, we show that besides the standard m-CAK 
critical point there exists a second one when the star's rotational speed is taken into account 
for a hot star with radiation driven wind. Using a simple /?— field, approximation we developed 
an analytical description for the location of the critical point and the value of the mass 
loss rate of the star. Our solution remains within a 1% confidence when compared with 
numerical calculations. This point is very important when trial solutions for numerical 
codes are selected. 

For the case of extreme high rotational speed (applied to a Be star), we have studied new 
solutions. Here exists only one critical point and numerical wind solutions were obtained. 
This new solution gives very slow terminal velocities (in the range ~ 190 kmsec^ 1 to ~ 450 
kmsec -1 , depending on the assumed values of the line force parameters) and the density 
of the wind is about 10 times denser than the standard m-CAK wind solution. These 
results indicate that new wind solution of the standard m-CAK model may be important 
in forming the disk in Be Stars. But to improve this model for Be stars, it is necessary that 
the following effects be to taken into account: the dependence of the radius and temperature 
as a function of rotational speed and latitude, a self-consistent treatment of the line force 
parameters, a detailed study of the influence of the boundary conditions, a re-analysis of 
angular momentum conservation (including viscous forces) and the (ignored here) influence 
of magnetic fields. 

Currently a detailed study of the topology of the singular points is underway. There are 
probably more solutions which involve more than one singular point as the case of magnetic 
fields. 

I would like to express my gratitude to N. Zamorano for his continuous support and 
critical comments on the manuscript. 

This work has been partially supported by Universidad de Valparaiso, internal project 
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DIUV 14/00 and 15/03. 



A. Coordinate Transformation. 

Here the basic steps toward equations (15a), (15b) and (15c) are outlined. The reader 
should keep in mind the original derivation by CAK. 



The partial derivatives of F(u,w,w/) (eq. 15a), with respect to u,w and wi are: 

l-h + <-- c '^ +CF ^)^ {mvirf (Al) 



9F A, , l\ .., „, fdCF , _CF g CF\ _, , , Q 



— = 1 + — U,' - C" -— + a 5 } g w~ d (ww/) a (A2) 

aw \ w 2 J \ ow w w J 

After using the new coordinate Y = w w' and Z = w/w', some derivative relation of 
the correction factor (see appendix B) and defining: 

dgiu) 

du 
where 



g(u) h(u) (A4) 



h{u) = S[-- . U ; I (A5) 



(A6) 



the singularity condition (w' dF/dw' = 0) reads: 

f 1 - JL) y - C f 2 (u, Z) g Z- & ' 2 Y^' 2 = 

here f 1 (u,Z) = CF(u,Z), 

f 2 (u, Z) = a /i(u, Z) - u Z e(u, Z) (A7) 

and e(u, Z) is a function defined in appendix B. 

The regularity condition (Z dF/du = 0) transform to: 



(A8) 



with 

/ 3 (u, Z) = (3m + Z) Z e(u, Z) + f\(u, Z) (h(u) Z + a - 5) (A9) 
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B. The Correction Factor 

The Correction Factor is defined by: 

2 f 1 ((l-i2 2 )v/r + n 2 v^ a 



CF 



i — A*2 



v' 



jidji. 



(Bl) 



where /j% — 1 — (R*/r) 2 and t>' = dv/dr. Integrating Bl and changing the variables from 
r,v to u,w (u — —R*/r, w = v/a, where a is the thermal velocity), the finite disc 
correction factor transform to: 



CF(u,w/w') 



1 1 



1 



l-au 2 (l 



1 w 



1 



W \ (1+a) 

1 — W — W 

W . 



:B2i 



where w' = dw /dr. Due to the fact, that CF depends on u and only the quotient Z = w/w', 
so defining A as: 

A = « (m + Z) (B3) 

we can re- write CF(u, Z = w/w') as: 



CF(A) 



1 1 



I- a A L 



1 - (1 - A) 



(l+a) 



(B4) 



The partial derivatives of CF with respect to u, w, w' are related to dCF/dX via the 

chain rule, namely: 

dCF dCF d\ 



du dX du ' 

dCF 8CF 8X 



dw dX dw 

dCF dCF dX 



dw' 



dX dw' 



Defining e(A) = e(u, A) = dCF/dX 



e(A) 



T-A) a -CF(A) 
X 



then (B5),(B6) and (B7) are related to (B8) by: 

1 dCF w'dCF 
e(A) 



w' dCF 



(B5) 
(B6) 
(B7) 



(B8) 



(B9) 



2u + Z du u dw uZ dw' 

Approximating Z by a /?— field (Z — (1 + u)/fl), we obtain CF and e only as functions of u. 
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Table 1. The location of the critical point, r crit , and the mass-loss rate, M, calculated by 
the analytical description for the CAK-model and compared against full numerical 

calculations. 
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Table 2. Analytical calculation for r cr u and M for the m-CAK model compared against 

full numerical calculations. 
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Table 3. Analytical calculation for r cr u and M for a CSPN compared against full 

numerical calculations. 
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Table 4. The value of the minimum rotational rate for the case when only one solution 

exists. 
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u=-R7r 



Fig. 1. — Behavior of the variable Z = w/w' as function of u (continuous line). Two /?— field 
approximations (straight-lines) for Z = (1 + u)/ (3 are over-plotted: j3 = 0.8 long-dashed 
line and f3 — 1.0 short-dashed line. 
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R(u,Z) 








Z=w/w' 



Fig. 2. — A view of R(u, Z) for typical values of u and Z for an Ob V star with v ro t/vb re akup 
0.5. 
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Fig. 3. — Same as fig. (2) but from a different perspective. 
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Fig. 4. — Same as fig. (3), but in addition of R(u,Z) the 'zero'-plane. The intersection 
of both surfaces gives the family of solutions that satisfies simultaneously the momentum 
equation and both singularity and regularity conditions. 
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R(u,Z)0 




Fig. 5. — Same as fig. (4), but here with v rot /vb re akup = 0.9. This figure shows clearly, that 
no solution exists for the standard m-CAK critical point for this rotational velocity. While 
a 'new' family of solutions is founded in the region where u < —0.2. 
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= r/R, 



Fig. 6. — Function f 123(11) vs. u for different values of (3 with line-force parameters: a = 
0.640 and 5 = 0.07. Continuous line: /3 — 1; dashed-dotted line: (3 = 2 and dashed line 
/3 = 3 
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Fig. 7. — Function Rcak(u) plotted against u for v rot /vbreakup = 0.0, continuous line; 0.1, 
dashed line; 0.3, dotted line; 0.5, dashed-dotted line; 0.7, long-dashed line. Figure a is for 
(3 = 1/2 and figure b is for /3 — 1. The location of the singular point is when Rcak{u) = 0. 
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0.92 



-0.9 



Fig. 8. — Behavior of Z = w/w' as function of u (continuous-line) in the interval [—1.0, —0.9]. 
The '+' symbol show the location of the critical point from full numerical calculations. Two 
/3-field approximations for Z — (1 +u)/(3 are over-plotted: (5 = 0.8 dashed-line and f3 = 1.0 
dotted-line. Clearly j3 = 1.0 intersects Ziu) almost at the critical point. 
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Fig. 9. — The function R app (u) vs. u for CAK (dashed-line) and m-CAK (continuous-line) 
models with (3 — 1 for different values of the rotational speed: a) v rot /vb re akup = 0.0; b) 0.5; 
c) 0.7; d) 0.9 for our 05 V test star. The roots of R app {u) gives the approximate location of 
the critical points. 
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Fig. 10. — In the phase-space u — Z, this figure shows a countourplot of R(u, Z) = 
(continuous-line) and a /3-field (dashed line, (5 — 1). The continuous-lines correspond to 
both critical point families, m-CAK and the new one. The parameters of R(u, Z) are the 
same as in figure (4). 
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Fig. 11. — The function R app {u) versus u for a Bl V star for two different values of the 
free-parameter (3; a) (5 — 1; b) (3 — 2.5. The different rotational speeds v rot /vbreakup are: 0.0, 
continuous line; 0.5, dashed line; 0.7, dotted line; 0.9, dashed-dotted line. 
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Fig. 12. — Velocity profiles v versus u for a B\ V star. Non-rotational case: dashed line; 
high rotational case with v rot /vbreakup = 0.8: continuous line. 
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Fig. 13.— Z versus u from full numerical calculations, continuous line; The 
approximation^ = (1 + u)//3 for (5 — 1 dashed-line, and for (3 = 2 dotted-line short-dashed 
line. 
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Fig. 14. — Function R(u, Z) from full numerical calculations, continuous line; R app (u) for 
(5 — 1 long-dashed line, and R app (u) for (3 = 2 short-dashed line. 
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Fig. 15. — Wind Density in g cm 3 versus u from full numerical calculations for the equatorial 
direction (continuous line) and polar direction (dashed-line). 



